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ABSTRACT 

Photo-heating associated with reionisation and kinetic feedback from core-collapse 
supernovae have previously been shown to suppress the high-rcdshift cosmic star for- 
mation rate. Here we investigate the interplay between photo-heating and supernova 
feedback using a set of cosmological, smoothed particle hydrodynamics simulations. 
We show that photo-heating and supernova feedback mutually amplify each other's 
ability to suppress the star formation rate. Our results demonstrate the importance 
of the simultaneous, non-independent inclusion of these two processes in models of 
galaxy formation to estimate the strength of the total negative feedback they exert. 
They may therefore be of particular relevance to semi-analytic models in which the 
effects of photo-heating and supernova feedback are implicitly assumed to act inde- 
pendently of each other. 

Key vifords: galaxies: high-redshift - galaxies: formation - galaxies: evolution 



1 INTRODUCTION 

The cosmic star formation rate (SFR) is an important ob- 
servable of our Universe. It is affected by a variety of phys- 
ical processes, many of which are in turn regulated by the 
SFR, giving rise to so-called feedback loops (for an overview 
see, e.g., Ciardi & Fcrrara 2005). Photo-ionisation heating 
due to the absorption of ionising photons from star-forming 
regions and the injection of kinetic energy from supernova 
(SN) explosions of massive stars provide two such feedback 
loops. Their implications for the assembly of the first gener- 
ation of galaxies have been extensively discussed in studies 
of t he epoch of reionisation (for a review of this epoch see, 
e.g.. lLoeb fc BarkanallioOll l. 

Photo-heating associated with reionisation increases 
the mean temperature of the interga l actic medium (IGM) 
to ~ 10* K (e.g., iHui fc Gnedh] 1 19971 ) a nd reduces 
the r ate at which hotter gas can cool (|Efstathioul 
1 19921 : IWiersma. Schave. fc SmithI l2009l ). The increase in 
the gas temperature keeps the IGM smooth and pre- 
vents the assembly of low-mass galaxies, that is, galax- 
ies with m asses corresponding to a yirial temperature < 
10* K (e.g., IShapiro. Giroux. fc Babullll994l : iGnedin fc Hul 
1 19981 ). Moreover, the gas in galcixies that have al- 
ready collapsed is relatively quickly photo-evaporated (e.g ., 
IShapiro. Iliev. fc Ragal l2004l : llliev. Shapiro, fc Ragal 120051 ). 
strongly decreasing the gas fraction of low-mass ha- 



los (e.g.. iThoul fc Weinberg! Il996l: iBarkana fc Loebl Il999l : 
iDiikstra et al.l 12004 ISusa fc Umemural I2006D . Indeed, the 
cosmic SFR has been predicted to exhibit a distinct dro p 
around the redshift of reionisation (|Barkana fc Loebl |2000| ). 
Photo-ionisation heating is therefore said to provide a neg- 
ative feedback on reionisationQ 

SN explosions of massive stars typically inject a few so- 
lar masses of gas with velocity of ~ 10* km s~^, correspond- 
ing to a kinetic energy of ~ 10^^ erg. The ejected material 
sweeps up and shock-heats the surrounding gas, entraining 
outflows sufficiently powerful to, at least temporarily, sub- 
stantially redu ce the gas fractions for galaxy-scale dark mat - 
ter halos (e.g.. lYepes et al.lll997l : IScannapieco et al.ll2006l ). 
Since this leads to a suppression of the SFR, SN explosions, 
like photo-heating from reionisation, provide a negative feed- 
back on reionisation. In addition to the depth of the gravi- 
tational potential, the ability of SN feedback to reduce the 
gas fractions gene rally depends on the geom etry of the gas 
distribution (e.g.. lMac Low fc Ferrarall 19991 ). 

Studies of the effects of photo-heating and SN feed- 
back on the SFR that considered each process in iso- 
lation have been augmented by studies that included 
both photo-heating and SN feedback. These studies in- 
clude simulation s of the forma-tion of the first stars (e.g . , 
iGreif et al.l [20071 : I Wise fc Abel l2008l : Iwhalen et al.ll2008i ). 
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^ As pointed out bv lPawlik. Schave. fc van Scherpenzeell | |2008|) . 
photo-heating also provides a strong positive feedback on reioni- 
sation because the increase in the gas temperature smoothes out 
density fluctuations, reducing the recombination rate. 
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of the evolution of isolated ga l axies (e.R., Fuiita et al.ll200-i : 



iDalla Vecchia fc Schayj |2008| . iTasker fc Bn^ 



2008fi 



and 



of ga laxies in a cosmological volume (e.g., iTassis et al.l 
I2OO3I ). Some of these studies also investigate the inter- 
play between_gh£to4ifia^ng_M SN feedback. For exam- 
ple, iKitavama fc Yoshidal ()2005l ) demonstrated in a one- 
dimensional hydrodynamical study that a previous episode 
of photo-heating may increase the efficiency of the evacua- 
tion of dark matter halos by (thermal) SN feedback and en- 
able the destruction of galaxies out to much larger masses. 
In this letter we report on another interaction of star for- 
mation feedbacks, using three-dimensional galaxy formation 
simulations. 

We employ a set of Smoothed Particle Hydrodynam- 
ics (SPH) cosmological simulations that include star forma- 
tion, and photo-ionisation heating from a uniform ultravi- 
olet (UV) background and/or kinetic feedback from core- 
collapse SNe. We investigate how photo-heating affects the 
high-redshift {z 6) SFR in the presence and absence of 
SN feedback and, conversely, how SN feedback affects the 
cosmic SFR in the presence and absence of a photo-ionising 
background. We find that the inclusion of SN feedback am- 
plifies the suppression of the cosmic SFR due to the inclusion 
of photo-heating. On the other hand, the inclusion of photo- 
heating amplifies the suppression of the cosmic SFR due to 
the inclusion of SN feedback. 

Photo-heating and SN feedback therefore mutually am- 
plify each other in suppressing the SFR. Our results are rel- 
evant to current implem entations of ( semi-) analytic models 
of galaxy formation (see lBaugh!l2006l for a review), in which 
the effects of photo-heating and SN feedback are implicitly 
assumed to act independently of each other. These models 
may thus underestimate the strength of the combined neg- 
ative feedback from photo-heating and SNe. 

We emphasise that none of our simulations has a suf- 
ficiently high resolution to achieve convergence in the cos- 
mic SFR. Future simulations will therefore be required to 
quantify our qualitative statement. We show, however, that 
the factor by which photo-heating and SN feedback amplify 
each other's ability to suppress the cosmic SFR becomes 
larger with increasing resolution, giving credibility to our 
main conclusion. 

This letter is organised as follows. We present our sim- 
ulation method in iJ2] and we illustrate our main result in 
^ Finally, we summarize our conclusions and discuss the 
caveats inherent to the present work in §3) 



2 SIMULATIONS 

Our simulation method is identica l to t hat employed in 
iPawlik. Schave. fc van Scherpenzeeil l|2008t ). to which we re- 
fer the reader for more details. We use a mo dified version o f 
the N-body/TreePM/SPH code GADGET-2 (|Springelll20"05l ') 
to perform a total of 12 cosmological SPH simulations at 
different resolutions, using different box sizes. We employ 
the set of cosmological parameters [57ni, f2b, f^A, eg, JT-s, /i] 
given by [0.258,0.0441,0.742,0.796,0.963,0 719], in agree- 
ment with the WMAP 5-year observations (iKomatsu et al.l 
|2008| ). The simulations include radiative cooling, star forma- 
tion and, optionally, photo-heating by a uniform UV back- 
ground and/or kinetic feedback from SNe (see Table [l]). 



Table 1. Simulation parameters. From left to right, table entries 
are: simulation label; comoving size of the simulation box, -Lbox! 
number of DM particles, A^dm! mass of dark matter particles, 
rrttjni . A prefix r9 indicates the inclusion of photo- heating in the 
optically thin limit from a uniform UV background (for ^ ^ 9) 
and a suffix winds indicates the inclusion of SN feedback (with 
initial wind velocity Hw = 600 km and mass loading 77 = 2). 
A bold font marks our set of reference simulations. 



simulation 


^box 




"T-dm 




[ h-^ Mpc] 




[10^ Mq] 


L6N256 


6.25 


256^ 


8.6 


r9L6N256 


6.25 


256^ 


8.6 


L6N256winds 


6.25 


256^ 


8.6 


r9L6N256winds 


6.25 


256^ 


8.6 


L3N256 


3.125 


256^ 


1.1 


r9L3N256 


3.125 


256^ 


1.1 


L3N256winds 


3.125 


256^ 


1.1 


r9L3N256winds 


3.125 


256^ 


1.1 


L3N128 


3.125 


128^ 


8.6 


r9L3N128 


3.125 


128^ 


8.6 


L3N128winds 


3.125 


128^ 


8.6 


r9L3N128winds 


3.125 


128^ 


8.6 



The gas is of primordial composition and is allowed 
to cool by coUisional ionisation and excitation, emission of 
free-free and recombination radiation and Compton cool- 
ing off the cosmic microwave background. Molecular hydro- 
gen is kept photo-dissociated at all times by the inclusion 
of a soft UV background. We employ the star formation 
recipe of ISchave fc Dalla Vecchial (|2008l '). to which we refer 
the reader for details. According to this recipe, gas forms 
stars at a pressure-dependent rate that reproduces the ob- 
served Kennicutt-Schmidt law ((Kennicutt 199 8). 

Photo-ionisation (heatin g) is included in t h e opt ically 
thin limit using a uniform iHaardt fc Madaul (|200lh UV 
background from quasars and galaxies for redshifts z ^ Zr = 
9. The value for Zr is consistent with the most recent deter- 
mination of the Thomson optical depth to wards reionisa- 
tion from the WMAP (5-year) experiment (|Komatsu et al.l 
l2008h . We inject an additional thermal energy of 2 eV 
per proton at z = z^. In the absence of shocks, gas par- 
ticles at the cosmic mean density are therefore kept at 
a temperature T » lO'^ K for z < z^ (see Fig. 1 of 
iPawlik. Schave. fc van Scherpenzeelll2008h . 

We model kinetic feedback from star formation us- 
ing the prescription of .Dalla Vecchia fc Schavd (|2008t ) , ac- 
cording to which core-collapse SNe locally inject kinetic 
energy and kick gas particles into winds. The feedback 
is specified by two parameters, the initial wind mass 
loading in units of the newly formed stellar mass, rj, 
and the initial wind velocity v^. We adopt i] — 2 and 
Vm = 600 km s~^, consistent with observat ions of local 
(e.g. IVeilleux" Cecil, fc Bland-Ha,wthornl 120051 ) and redshift 
z » 3 fe.g. IShaplev et al.l |2003| ) starburst galaxies. This 



choice of parameters implies that 40 per cent of the energy 
available from core-collapse SNe is injected as kinetic en- 
ergy (assuming a Chabrier initial mass function in the range 
0.1 s; M/Mq < 100 and 10^^ erg per M > 6 Mq star), while 
the remaining 60 per cent are implicitly assumed to be lost 
radiatively. 
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Figure 1. Top panel: evolution of the total SFRs in our set of 
reference simulations. Photo-heating and SN feedback were in- 
cluded for the red and dashed curves, respectively. Star forma- 
tion is strongly suppressed when photo- heating and/or kinetic 
feedback from SNe are included. Bottom panel: Suppression fac- 
tors and Sw^h (right y-axis) and feedback amplification factor 
X = Swh/C^wSh) (left y-axis). The fact that X > 1 implies that 
photo-heating and SN feedback amplify each other's ability to 
suppress the SFR. 



We use a Friends-of-Friends halo finder with linking 
length b = 0.2 to obtain a list of dark matter hales con- 
tained in each of our simulation outputs. Only halos con- 
taining more than 100 dark matter particles are included 
in these lists. For each simulation output we compute the 
SFR associated with dark matter halo as follows. First, gas 
particles are attached to the nearest dark matter particle. 
Second, the SFR of a dark matter halo is the sum of the 
SFRs of the gas particles that were attached to the dark 
matter particles it contains. 



3 RESULTS 

Our main result is shown in Fig. [T] The top panel shows 
the evolution of the total SFR, defined as the sum of 
the star formation rates over all halos, for our set of ref- 
erence simulations. We use p*,wh, p*,h and p,,w to de- 
note, respectively, the SFR densities in the simulations with 
both SN feedback and photo-heating (red dashed curve, 
r9L6N256winds) , with photo-heating but without SN feed- 
back (red solid curve, r9L6N256) and with SN feedback but 
without photo- heating (blue dashed curve, L6N256winds). 
We denote the SFR density in the simulation that included 
neither photo-heating nor wind feedback with p, (blue solid 
curve, L6N256). 

Both the inclusion of photo-heating and the inclusion of 
kinetic feedback from SNe lead to a significant suppression 
of the SFR. The factor Sh = p./p,,h by which the SFR is 
suppressed due to photo-heating is smaller than the factor 
Sw = P*/p*,w by which the SFR is suppressed due to kinetic 
feedback. As expected, the simultaneous inclusion of photo- 
heating and feedback from SNe leads to a suppression of 
the SFR by a factor Swh = P*/p*, wh that is larger than the 
factors Sh and Sw by which the SFRs are suppressed due to 
the sole inclusion of either photo-heating or SN feedback. 



Interestingly, the factor by which the SFR is suppressed 
due to photo-heating is larger in the presence (set of dashed 
curves) than in the absence (set of solid curves) of SN feed- 
back. Conversely, the factor by which the SFR is suppressed 
due to SN feedback is larger in the presence (set of red 
curves) than in the absence (set of blue curves) of photo- 
heating. 

Photo-heating and SN feedback thus mutually am- 
plify each other in suppressing the SFR. This amplifi- 
cation probably arises because the inclusion of photo- 
heating keeps the gas diffuse, which makes it easier for 
the winds to drag it out of halos. Winds, on the other 
hand, move gas from the central to the outer parts of ha- 
los, where it is more susceptible to the photo-evaporation 
process. Models that implicitly ignore this interaction be- 
tween photo-heating and SN feedback, like for exam- 
ple (semi-) analytic mode ls of galaxy formation (e.g., 
Khochfar fc Ostriked boOSi: [Monaco. Fontanot. fc Taffonil 



2007 ; iBenson et aLl liooa ICroton et al.l l2006l : isomerville 
20021 ). thus may underestimate the strength of the feedback 
these processes exert. 

We now define the feedback amplification factor x = 
Swh/(swSh). A value X = 1 would indicate that photo- 
heating and SN feedback suppress the SFR independently of 
each other. A value X > f (x < 1) would indicate that photo- 
heating and SN feedback amplify (weaken) each other's abil- 
ity to suppress the SFR. The evolution of the amplification 
factor X is shown in the bottom panel of Fig. [T] together 
with that of Swh and s„Sh. The fact that x > 1 for z < 9 
implies that photo-heating and SN feedback amplify each 
other in suppressing star formation^ 

To demonstrate that this amplification is indeed mu- 
tual, we write s^h = s^ihSh = s^^s„. This defines the sup- 
pression factors s„|h = p*,h/p.,wh and Sh|w = P«,w/p*,wh. 
Thus, Sw|h is the factor by which SN feedback suppresses 
the SFR in simulations that include photo-heating and s^iw 
is the factor by which photo-heating suppresses the SFR in 
simulations that include SN feedback. We have 

P»,h/p*,wh _ P»,w/p *,wh _ £hh 

Sh 



Sw|h 

Sw 



p./p., 



p*/p* 



(1) 



This shows explicitly that the amplification of the suppres- 
sion of the SFR due to photo-heating in simulations that 
include SN feedback is equal to the amplification of the sup- 
pression of the SFR due to SN feedback in simulations that 
include photo-heating. It implies that we cannot determine 
whether photo-heating amplifies the effect of SN feedback 
or vice versa. 

Fig. [2] shows, for our set of reference simulations, the 
dependence of the suppression of the SFR due to photo- 
heating and/or SN feedback on halo mass. The top panel 
shows the cumulative SFR at z = 6 in dark matter halos 
less massive than Mdm- The bottom panel shows, similar to 
the bottom panel of Fig.[T] the feedback amplification factor 
X obtained from analogous definitions of the suppression fac- 
tors Sh, Sw and s„h applied to the cumulative SFR at z = 6 
shown in the top panel. The vertical dotted lines indicate 
the dark matter mass corresponding to a virial temperature 
Tvi, = 10* K. 

^ Note that in our simulations Sh = 1 = X foi' 2^9 because 
there is no photo-heating at these redshifts. 
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Figure 2. Top panel: cumulative SFRs as function of dark matter 
halo mass at 2 = 6 for our set of reference simulations. Photo- 
heating and SN feedback suppress the SFR over roughly comple- 
mentary mass ranges. Bottom panel: similar to the bottom panel 
of Fig. [T] but now for the cumulative SFRs shown in the top 
panel. The vertical dotted lines indicate the dark matter mass 
corresponding to a virial temperature Tvir = 10** K. 



Figure 3. Evolution (left-hand panel) and mass-dependence (at 
2 = 6; right-hand panel) of the feedback amplification factor x 
for different choices of the box size and resolution. The vertical, 
dotted line in the right panel indicates the dark matter mass 
corresponding to a virial temperature Tvir = 10* K. The feedback 
amplification factor is insensitive to the size of the simulation box 
but increases strongly with resolution. 



While photo-heating strongly decreases the SFR in ha- 
los with virial temperatures Tvir < 10* K, the inclusion of SN 
feedback leads to a strong decrease in the SFR in halos with 
Tvir 10* K. Radiative and kinetic feedback thus act mostly 
over complementary mass ranges. This dichotomy arises be- 
cause photo-evaporation mainly affects the gas fractions of 
halos with masses that correspond to virial temperatures 
that are of order of or smaller than the thermal temper- 
ature to which the gas is photo-heated, Tvir < 10* K. In 
contrast, star formation and the associated SN feedback be- 
come only efficient for halos with masses corresponding to 
Tvir ^ 10* K, because our simulations do not include ra- 
diative cooling from metals and molecules (and because of 
our limited resolution), which could potentially enable star 
formation in halos of much smaller virial temperatures. 

The left-hand (right-hand) panels of Fig. |3] show the 
evolution (mass-dependence) of the feedback amplification 
factor X obtained from simulations for which we have varied 
the size of the simulation box and/or the resolution. Note 
that the solid, black curves in the left-hand and right-hand 
panels are identical to the solid, black curves shown in the 
bottom panels of Figs. [T] and (2] respectively. 

Changing the size of the simulation box by a factor 
of two (at fixed resolution; solid curves) has little effect. 
On the other hand, increasing the resolution by a factor of 
two (while keeping the size of the simulation box fixed; red 
curves) significantly increases x- The increase in x with in- 
creasing resolution is likely due to both an increase in the 
fraction of galaxies that are subject to photo-evaporation 
and an increase in the SFR (and thus associated SN feed- 
back) of all galaxies. 



4 DISCUSSION 

Photo-heating from reionisation and supernova (SN) feed- 
back are key processes that determine the star formation 



rate (SFR) in the high-redshift Universe. Using a set of 
cosmological SPH simulations that include radiative cooling 
and star formation, we analysed the z > 6 star formation 
history in the presence of photo-ionisation heating from a 
uniform UV background and/or kinetic feedback from core- 
collapse SNe. The inclusion of photo-heating and SN feed- 
back both lead to a suppression of the SFR. 

We showed that the factor by which the SFR is sup- 
pressed due to photo-heating is larger in the presence than 
in the absence of SN feedback. We also showed that the fac- 
tor by which the SFR is suppressed due to the inclusion of 
SN feedback is larger in the presence than in the absence 
of photo-heating. This mutual amplification of SN feedback 
and reionisation heating is the central result of the present 
work. 

We caution the reader that our simulations have not 
fully converged with respect to resolution and that we have 
ignored some potentially important physical processes. This 
result will therefore need to be confirmed and quantified 
with future simulations. In what follows we briefly discuss 
the most important physical effects that our analysis ig- 
nored. 

We computed photo-heating rates from a uniform UV 
background in the optically thin limit. Our simulations 
therefore do not account for the self-shielding of gas 
from ionising radiation. This may lower the fraction of 
the ga s that is photo-e va porat ed (e.g.. [ Kitavama et al.l 



2000l: ISusa fc Umemural (20041 : iDiikstra et all |2004| ) 



Iliev Shapiro, fc Ragal l2005l ) (extending the work of 



Shapiro. Iliev. fc Ragal |2004| ) have, however, pointed out 



that regions that are initially self-shielded will eventually 
also be photo-evaporated: as subsequent layers of gas are 
photo-evaporated, previously self-shielded regions become 
exposed to ionising radiation and are eventually stripped 
away. Moreover, the evaporation of the initially self-shielded 
gas may proceed at a speed comparable to the speed pre- 
dicted by simulations that compute photo-heating rates in 
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the o ptically thin limit (Fig. 3 in llliev. Shapiro, fc Ragal 
I2OO5I ). Note also that because SN explosions decrease 
the gas density, the effects of self-shielding may be less 
prominent in simulations that include this type of feed- 
back. Clearly, the importance of self-shielding and the 
applicability of the optically thin limit to the present 
problem need to be critically assessed using cosmological 
radiation- hydrodynamical simulations. 

Although all of our simulations employ a sufficiently 
high resolution to resolve all halos with virial temperatures 
Tvir ^ 10* K with at least 100 particles, none of them has the 
resolution to properly reproduce the properties of the multi- 
phase medium associated with the star-forming regions these 
halos host. The factors by which photo-heating and SN feed- 
back suppress the SFR are therefore not yet converged. We 
have, however, demonstrated that the effect of mutual am- 
plification of photo-heating and SN feedback becomes only 
stronger with increasing resolution. 

Because we have assumed the presence of a photo- 
dissociating background, the formation of molecular hydro- 
gen is suppressed in our simulations. In reality, the gas 
may contain a significant frac tion of molecular hydrogen be - 
fore reionisation (but see, e.g.. lHaiman. Rees. fc Loeblll997l ). 
Star formation and the associated kinetic feedback would 
then already be efficient in halo s with virial tempera tures 
much smaller than 10* K (e.g.. iTegmark et al.l riggTl ). Be- 
cause the gas fraction in these halos would be affected by 
both photo-heating and SN feedback, we expect that the 
inclusion of molecular hydrogen would only strengthen our 
main result. 

We have also ignored the existence of atoms and ions 
heavier than helium. SN explosions may, however, quickly 
enrich the interstellar and intergalact ic gas with metals (e.g., 
iBromm. Yoshida. fc Hern^^l2003l ). which would increase 
its ability to cool. Metal enrichment thus provides a posi- 
tive feedback that may partially offset the negative kinetic 
feedback from SN explosions. 

None of the caveats we discussed seems, however, likely 
to invalidate our main qualitative conclusion that photo- 
heating and SN explosions amplify each other's effect on 
the cosmic SFR. Galaxy formation models that treat the ef- 
fects of photo-heating and SN explosions independently of 
each other, like e.g. (semi-) analytic models, may therefore 
significantly underestimate the effect of these feedback pro- 
cesses on the SFR. Because photo-heating and SN feedback 
are important processes that affect the gas in low mass-halos 
at all epochs, our findings may also have implications for the 
understanding of the properties of low-redshift galaxies, e.g. 
in the context of th e missing satellite pr oblem (for a recent 
discussion see, e.g., IKoposov et aLll2009[ ). 
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